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Abstract 

We consider the non uniform AdS black string equations in arbitrary number of dimension in a per- 
turbative approach up to order 2 and in a non perturbative. We restrict the study in the perturbative 
approach to the backreacting modes, since they provide the first relevant corrections on the thermo- 
dynamical quantities of the solutions. We also present some preliminary results in the construction of 
non-perturbative solutions, in particular, we present a first part of the non uniform - uniform black 
string phase diagram. Our results suggests the existence of a new stable phase for AdS non uniform 
black strings, namely long non uniform black string, with the extra direction length of the order of the 
AdS curvature. 

Keywords: Black strings, non uniform, AdS/CFT, higher dimensional gravity. 

1 Introduction 

There has been a growing interest for solutions in AdS spacetime during the last decade, largely motivated 
by the AdS/CFT correspondence [2 Hj. This correspondence relates the geometry of an asymptotically 
AdS space with a conformal field theory defined on the background of the conformal boundary of the AdS 
space. In this context, black objects are of great interest; for example, the Hawking-Page phase transition 
[3] between the five dimensional spherically symmetric black holes and the thermal AdS background was 
interpreted by Witten, through AdS/CFT, as a thermal phase transition from a confining to a deconfining 
phase in the dual four dimensional M = 4 super Yang-Mills theory [3] . 

In four dimensional asymptotically flat space, the uniqueness theorem guaranties that the horizon topol- 
ogy of black objects is 52- In higher dimensions, this theorem does not hold and several black objects with 
various topology have already been constructed. Among them, let us mention the black strings [5] of horizon 
topology Sds x Si, black rings [5] (toroidal horizon topology), black holes (Sd-3 horizon topology) [7] and 
their generalisation with charge, rotation, cosmological constant, etc. 

In 1993, Gregory and Laflamme [5] showed that black strings are unstable toward long wavelength 
perturbations in asymptotically locally flat space (see [9] for a review). Just after that discovery, it was 
widely believed that the end point of the black string instability in asymptotically locally flat spacetime 
should be a caged black hole, where the asymptotical space in d + 1 dimensions is Aid x Si (with Aid the 
e?-dimcnsional Minkowski space), but it was argued [10| that the transition between a black string and a 
caged black hole would take an infinite proper time at the horizon and thus even 'more' for an observer 
lying outside the horizon. This motivated the study of extended black objects which arc not translationally 
invariant along the Si coordinate, namely the non uniform black strings, first in a perturbative approach 
[TT] (see also [T^]) and later in the full non linear regime [T3]. The resulting thcrmodynamical phase diagram 
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|14j is now well known although dynamical phase transitions are still to be constructed. In this diagram, the 
non uniform and the caged black hole branch meet at the merger point, where a topological phase transition 
is expected to occur. 

More recently, the question of the black strings instability in asymptotically locally AdS spacctime has 
been addressed and it has been argued that this instability should persist in asymptotically locally AdS space 
in any number of dimensions for small AdS black strings [15| where the ratio r/ l /£ << 1, with the horizon 
radius and £ the AdS radius, fixed by the cosmological constant A. Note that there exists another phase of 
AdS black strings, namely large AdS black strings [H1[T7], where r^/i > 1 which are thermodynamically 
and dynamically stable |15] . As already noticed in [3] , AdS space acts like a confining box: when the ratio 
rh/£ becomes larger than some critical values, the wavelength of the instable modes cannot fit the 'AdS box' 
and thus the configuration is stable. 

However, in asymptotically locally AdS spaces, much less is known about the counterpart of the solutions 
of A = Kaluza-Klcin black objects. The evolution of AdS black holes has been studied in [TB] but still has 
to be done for AdS black strings. The thermodynamical properties of the uniform phase is now well known 
|16j and non uniform solution have been constructed in |19j in a perturbativc approach and in 6 dimensions. 
In order to construct a phase diagram in AdS, one has to consider various stationary solutions and compare 
their thermodynamical properties since these stationary solutions should be the equilibrium configurations. 
This could give some clues on the endpoint of dynamical evolution and phase transitions between different 
black objects in AdS. 

In this paper, we consider perturbative non uniform black strings in asymptotically locally AdS spacetime 
for arbitrary number of dimensions. The decay of the non uniform modes in arbitrary number of dimension 
is sharp so it is difficult to construct the global conserved quantities. Instead, we focus on the horizon 
thermodynamical properties such as the entropy and the Hawking temperature in order to investigate the 
possible existence of new stable phases of non uniform black strings; namely long black strings, which should 
be counterparts of big AdS black strings. It will turn out that this stable phase indeed exists and appears 
only in the non uniform phase. We also construct the non perturbativc non uniform black string solution 
for some different values of the cosmological constant. 

This paper is organised as follows: in section 2, we present the model. In section 3, we give the near 
horizon and asymptotical solution to the equations while we present the boundary conditions and numerical 
technique in section 4. Section 5 is devoted to the thermodynamical properties of the solutions, we present 
the numerical results and give a discussion of the solutions in section 6. Finally we present the solution to 
the non perturbative problem in section 7 and discuss the properties of the non uniform black string in the 
regime where the cosmological constant is small. The results are summarized in the conclusion. 



2 The model 

We consider the c?-dimcnsional Einstcin-Hilbert action supplemented with the Hawking-Gibbons boundary 
terms and a negative cosmological constant: 

where G is the d-dimcnsional Newton constant, g is the determinant of the metric on the spacetme manifold 
M., h the determinant of the induced metric on the boundary manifold dM., R is the scalar curvature, £ is 
the AdS radius and K is the trace of the extrinsic curvature of the boundary spacetime. In what follows, 
we will work in geometrical units where G = 1. 

We supplement the action (TTJ) with the following ansatz, which is relevant in the study of non uniform 
black strings [11] : 

ds 2 = -b(r)e 2A ^dt 2 + e 2B ^ (J^- + a(r)dz 2 ^ + e 2C ^r 2 dQ 2 d _ 3 , (2) 
where df2^_ 3 is the square of the line element on the unit (d— 3)-sphere and z£ [0, L], L being a real number. 
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The ansatz ^ leads to partial differential equations which are technically difficult to deal with. In order 
to obtain a set of ordinary differential equations, we develop the non uniformity (the z depending part) in 
a Fourier series and in term of a small parameter e: 

X(r, z) = eXi(r) cos(fcz) + e 2 (X (r) + X 2 (r) cos(2kz)) + 0(e) 3 , (3) 

where X generically denotes A, B, C and k = 2tt/L. 

With the expansion ([3]), we obtain ordinary differential equations at each order in e and for each indepen- 
dent Fourier mode. The equations at order e° are the background equations; the solutions of these equations 
have been studied in [IB]. The order e 1 leads to the linear stability of the background and the solutions have 
been constructed in [T5] , 

In order e 2 , there are two independent modes: the Xq and the X 2 fields. The Xq modes are the 
backreacting fields and have been studied in [T9j for d — 6. They arc relevant for the first corrections on 
thcrmodynamical quantities of non uniform black strings. The X 2 modes are massive modes and do not 
contribute in thermodynamical corrections (at first relevant corrections In addition, the X 2 modes are 
independent of the Xo fields since they are different Fourier modes; that's the reason why we will consider 
only the X fields at order e 2 in this paper. 

The equations for the background can be found in [16] and the equations for the first order are in |15j . 
The equations for the second order in e are quite long and we prefer not writing them here but they are 
straightforward to obtain. Let us just remind that it is possible to solve algebraically for the first order 
field B\ (resp. second order field Bo ) as a function of the other first order field (resp. second order fields). 
Let mention also that the equations for the backreacting fields present the same shift invariance than in the 
6-dimensional case, Xq — > Xq + const. 



3 Near-horizon and asymptotic expansion 

In order to compute the asymptotic expansion, two cases must be distinguished: the case of odd number 
of dimensions and the case of even number of dimensions. In odd number of dimension, log terms arise in 
the background asymptotic expansion while in even number of dimensions, there are no fog terms in the 
expansion; these terms might have a nontrivial repercussion on the first and second order fields, but it turns 
out that at the leading order of the asymptotic expansion, it is not the case. 

The background fields obey the Fcffcrman- Graham expansion which can be found in |16j for the case of 
interest. Let us just give the leading order expansion of the background fields: 

o(r)w6(r)«/(r) = J + (4) 

which is the appropriate asymptotic behaviour of a locally asymptotically AdS spacctime. 

In such an asymptotic background, the leading order of the asymptotic expansion for the first order fields 
is given by: 



A x {r) = _(d-3)7i^-J +0 {-) > C i( r )=K[-) + °{r) ' (5) 

with 71 a real constant to be determined numerically. The influence of the log terms we mentioned shows up 
in higher order terms in the expansion and is not relevant for our analysis. In fact, this is not the most general 
asymptotic expansion for the first order fields. The most general expansion contains constant terms and 
terms of order r~^ d ~ 4 \ These terms cancel once the boundary condition is imposed (Ai(r) — > 0, C\{r) — > 
when r — > oo; see next section). 

The backreacting fields A , Co follow the same pattern: 

/ p\ d ~ 1 / f\ d + 2 /f\ d ~ 1 / f\ d+2 

Ao(r) = -(d-3)7o(-J + {~) > 0)(r)=7o(;J +Of-J , (6) 



1 thermodynamical quantities involve an integration over z from to 2tt/L, thus suppressing X2 terms, see 1191 . 
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where 70 is a real constant, also to be determined numerically. Once again, in the most general expansion, 
there are lower order terms which would give infinite contribution to the mass and thus are unphysical |19j . 

Higher order terms for first and second order fields can be obtained after a straightforward calculation. 

The near horizon behaviour of the background fields and first order are reminded here: 



2a h (d-l)r h 2a h (d-l) 2 rl (r - r,) 2 

air) = ah + t~, tt^ n — 7T-?{ r ~ r h) H o 7: \-U[r — rh) 

(d-A)P + (d-l)rl y ((d-i)P + {d-l)r}f 
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b h ( d-4)((d-3)£ 2 + (d-l)r 2 ) {r-r_ h f 
(d-4)£ 2 r h + (d-l)r 3 h 2 



6(r) = b h (r-r h )+ -^~,/ , ^ ll2 _"' i ' f J ^J'' n,v ^' + 0{r-r h f, (7) 

rf 

Ai(r) = A 1Q + A n (r - r h ) + 0(r - r h ) 2 , d(r) = C 10 + C n (r - r h ) + 0(r - r h ) 2 , (8) 
where 



_ | (^ + (^\ (r _ ri) _ ( ,_ 4) ^ + ^(^ i )! +0(r _ rk) . 



An = 



2a, (Am - do) (d - 4) (d - 3) ^ 2 + (2a, A 10 (d - 5) (d - 1) + (-2A W + C 10 (d - 3)) k 2 £ 2 ) r 2 



h 



3a h r h ((d-4)£ 2 + (d-l)r 2 ) 



2 (Aw ~ C10) C10 (2a,(d-l) + fc 2 l 2 )r, 
11 r, a,(d-4K 2 + a,(d-l)r 2 ' W 

and a,, 6, are normalisation constants to be fixed such that the background fields follow j4 10 ,Cio arc 
real constants. Note that due to the linearity of the first order equation, either A w or C\q can be fixed 
arbitrarily, only the ratio A\q / C10 is non arbitrary. 

The second order fields have the following near-horizon expansion: 

Mr) = A aa + A 01 (r - r h ) + 0(r - r h ) 2 , C (r) = C 00 + C i (r - r h ) + 0(r - r h f, (10) 

where Aoo, Coo are to be fixed using the invariance under Aq — > Aq + const, Co — > Co + const such that the 
fields Aq, Co decay to at infinity, Cn is an arbitrary real constant and 

a h C Q1 (d - 4) (d - 3) I 2 + (A 2 10 + 2A 10 C 10 (d - 3) - C 2 (d - 3)) k 2 £ 2 r h + a h C 01 (d - 3) (d - 1) r 2 
01 ~ 3a, (d - 4) £ 2 + 3a, (d - 1) r 2 

This expression for Aqi is obtained by imposing regularity of the equation in the near horizon region; if 
Aqi is not chosen to be the above expression, there will exist a diverging term proportional to (r — in 
the near horizon limit. 



4 Boundary conditions and numerical technique 

The numerical technique follows that of reference |19j : we first integrate the background fields with appro- 
priate boundary conditions (see ref. |16| ) then we integrate the first order fields (see |15j for the boundary 
conditions). The backreacting fields are then integrated with the following initial conditions: 

A (r h ) = a , C(r h ) = X o, C'(r h )= X i (12) 
OhXi (d - 4) (d - 3) £ 2 + (A 2 W + 2A 10 C 10 (d - 3) - C 2 (d - 3)) k 2 £ 2 r h + a hXl (d - 3) (d - 1) r 2 
° [rh) 3a,(d-4)£ 2 +3a,(d-l)r 2 

where ao, Xo ar c arbitrary constant to be fixed a posteriori using the shift invariance A)(r) — > Ao(r) + const, 
Co(r) — > Co(r) + const, while xi 1S a constant which is tuned such that the fields Aq,Cq follow the decay 
©■ 

In practice, we integrate the background and the first order with the solver Colsys |20j and the second 
order is integrated using a Runge Kutta algorithm at order 4. The integration is carried out from the horizon 
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radius, ru to some R >> r/j. The fields Aq, Co follow the asymptotic © if RC' (R) + (d — 1)Cq(R) = 0. 
The problem here is that Co can always be shifted by an arbitrary constant so we impose the decay (jSJ) 
transposed to the derivative of Co: 

RC \R)+dC (R) = 0, (13) 

the value of C' '(R) being obtained using the field equations. 

The value of RCq(R) + dC' (R) is a function of Xij say C(xi)- The backreacting fields will follow the 
asymptotic decay ^ if xi is chosen to be a root of C(xi)- We used the Newton algorithm in order to find 
the value of xi such that C(xi) = 0. 



5 Thermodynamical properties 

The entropy is given by a quarter of the event horizon area and the first relevant correction to the entropy 
coming from the non uniformity appears at the backreacting level. Let us write the entropy as S = Sq + £ 2 5S, 
So being the entropy of the uniform black string. 8S the correction from the backreaction. Then, 

So T i (14) 

- ota a h (-Ala + 2A 10 C 1Q (d - 3) + (4C 00 + C? (d - 3)) (d - 3)) (d - 1) + C 10 (A 10 + C w ){d - 3)k 2 £ 2 
SS/S ° ~ 4a h (d-l) ' 

where is the surface of the unit (d — 3)-sphere and the background quantities are evaluated with the 

length L = 2n/k, k being determined from the solution to the first order equations. 

The temperature at the horizon can be computed demanding regularity in the Euclidean sections and is 
given by: 



T * = iv^ ((d_lH+(d_4)n (15) 

2a (d - 1) (2A,o + A\ Q ) + (d - 3)(A 10 - C w )C w k 2 £ 2 
6Th/TH = 4^1) ' 

where Th is the background temperature and STh the corrections from the backreaction, the temperature 
of the non uniform phase being denoted by Th + £ 2 8Th- 

The mass and tension can be computed using the counterterm procedure j'21j . We refer the reader to 
[TBI rjS] for the details of the technique and an application to uniform asymptotically AdS black string 
and perturbative non uniform AdS black strings in 6 dimensions. Let us notice that there exists another 
procedure, free of ambiguities that might arise during the standard counterterm procedure, namely the 
Kounterterm procedure [22]. 

However, the thermodynamical quantities involved in these procedures are difficult to extract from the 
numerical solution in the general d-dimensional case because of the backreacting fields asymptotical decay 
([5]) . Instead, it is in principle possible to compute the mass by integrating the first law with fixed asymptotical 
length (L fixed) and to use the Smarr formula [16j to extract the tension. In practice, this procedure works 
well except for very small value of the horizon radius (with AdS length fixed to one), which is the region of 
interest in the perturbative approach. We shall come back to this point in the discussion. For our purpose, 
it is sufficient to evaluate only the entropy and Hawking temperature (with corrections arising form the 
backreactions) . 

In order to investigate the thermodynamical stability of the non uniform phase, we consider the specific 
heat, defined as the derivative of the entropy with respect to the temperature. In the case of perturbative 
non uniform black strings at order e 2 , the specific heat reduces to 
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with SS and STh as defined above. 

This quantity gives the variation of entropy whith respect to the Hawking temperature in the non 
uniform phase at the merger point; non uniform solutions with C p < are thermodynamically unstable, 
while solutions with C p > arc thermodynamically stable. Consequently, in order to investigate the existence 
of a stable phase, it is sufficient to compare the sign of SS and STh- 

6 Numerical results and discussion 

We integrated the field equations for values of d from 5 to 15. We considered only perturbative non uniform 
black strings with critical length. The critical length is well defined only for small AdS black strings (for big 
AdS black strings, the critical wavenumer k becomes imaginary and so does the critical length L = 2n/k). 
Our results show that for every number of dimensions considered, there is a phase of non uniform black 
string presenting a positive specific heat (Fig. [T]). 




Figure 1: The correction on the entropy SS as a function of the correction on the Hawking temperature STh 
for various number of dimensions. There arc points on this figure with STh > corresponding to large value 
of This stable phase occurs for small value of k (i.e. large value of L) and is the analogue of the big 
AdS black string, with the large length being now in the S 1 direction. The small figure in upper left box is 
a zoom of the region with STh > 0. 



This phase arises for small values of the critical wave number, i.e. large value of the length (with £ fixed). 
This is the analogue of the small-big black strings or black hole in AdS where the relevant parameter for 
the thcrmodynamical stability was r^jl. Here, the relevant parameter is [i\ — L/£; short non uniform black 
strings solutions with fix « 1 arc thermodynamically unstable while long non uniform black strings with 
fj,i « 1 are thermodynamically stable. 

Note that this effect is not present for the uniform phase. In the uniform phase, the length enters only 
as an overall factor in the thermodynamical quantities. In the non uniform phase, the thermodynamical 
quantities depend on the length in a non trivial way (via k = 2t:/L). Let us stress the fact that these 
short-long black strings are present in the small AdS phase, since it is the phase we are dealing with. We 
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Tabic 1: The value of the ratio ix\ s = L/£ where the stable phase occurs for various d. 



expect this property to be a generic feature of non uniform black strings in arbitrary number of dimensions. 

Figure [2] shows the direction of the phase transition in a S — Th diagram for d = 9. In the short black 
string phase (i.e. with L << £), the entropy increases and the temperature decreases along the non uniform 
phase while in the long black string phase (L £) the entropy increases and so does the temperature. Let 
us emphasize once again that this new phase of non uniform black strings occurs for small black strings 
(r « £) since the length is not well defined for big black strings (r « £). 




0.7 0.8 0.9 1.0 1.1 1.2 



Figure 2: The phase diagram in the S — Th plane for d = 9. The background and the direction of the 
corrections for some value of the rescaled length are shown. The same pattern holds for all the number of 
dimension we considered. 



We wanted to compute the mass and tension for the solutions we constructed in order to investigate the 
existence of a (^-dependinjQ) critical dimension, but it turned out that these quantities are much difficult 
to evaluate precisely. The problem comes from the small values of the horizon radius: when integrating 
the background, we impose generically a(r^) = 1, b(rh) = f(rh) = 0, and use the arbitrariness of the a, b 
functions normalisation in order to impose a, b, f — > r 2 /£ 2 asymptotically, a posteriori. For small value of r^, 
the asymptotical region appears for larger and larger value of the radial coordinate; this introduces numerical 
noise in the determination of the normalisation factors for the background functions and thus affects the 
horizon quantities. Moreover, the small black string phase is defined for < 1 (we are working with £ = 1; 
fixing another value of I doesn't solve the problem.). The investigation of a critical dimension thus deserves 
further attention. 

However, we expect this long non uniform black string phase to be relevant in the potential existence of 

2 with fixed horizon radius or equivalently rrj depending, with fixed I. 
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a ^-depending critical dimension, since there is a major change in the thermodynamical properties of this 
phase. 



7 Non Uniform Black String 

In this section we present some preliminary results on the solution to the full non-linear system of partial 
differential equations in six dimensions. These equations can obtained starting from the ansatz without 
developing the non uniformity as a Fourier Series. However, demanding regularity on the horizon leads to 
non trivial relations between the functions A, B, C. 

Changing the parametrisation of the background (where A = B = C = 0)to& conformal-likc gauge, 

ds 2 = -b{f)dt 2 + + ~a{f)dz 2 + (r 2 + f 2 )d^ (17) 

(r 2 h + r 2 )f(r) 

and using a non uniform ansatz in this gauge, 

ds 2 = -b(r)e 2 ^dt 2 + e 2 ^ f / 2 f; 2 - + d(r)dz 2 ) + (r 2 + f 2 )e 2 ^^ 2 _ 3 , (18) 

leads to much simpler regularity conditions at the horizon, which will be given later. 

In this parametrisation, the horizon is located at f = and the functions b and / behave like f 2 close to 
the horizon. Note that the variable f and the functions 5, b, f are related to the variable r and the functions 
a,b, f according to 



r = \Jrl + f 2 , 5(f) = a(r) , b(f) = b(r) , /(f) = f(r). (19) 
The equations with the parametrisation (|18p are given by 

_ 5e 2B(f,z) r 2 £'(f) a'(f)b'{f) b'{ff 6'(f)/'(f) b'(f)g'{f) 



s 2 /(f) 5 (f) 2rb(r) 4a(f)b(f) 46(f) 46(f)/(f) b(f)g(f) 

b"(f) r 2i(oa)(^ z f 3r 2 i(°- 1 )(f,z)C^ 1 )(f, z) r 2 A^(r,z) A^°\f,z) 
+ 26(f) + L 2 5(f)/(f) 5 (f) + L 2 a(r)f(r)g(r) + L 2 5(f)/(f) ff (f) r 

a'{f)A^{f,z) b'{r)A^(r,z) f'(f)A^°\f,z) 2 5 '(f)i( 1 -°'(f, z) 
+ 25(f) + b(f) + 2/(f) + W) 

+ A™ (f , z f + g^glg^fe) + 34(1,0) (r ~ ; ^£(1,0) ( - z) + 4(2,0) ( - ) z) = Q) 
26(f) 



3e 2B(f, Z )-2C(f,») r 2 5e 2B(f, Z ) r 2 ^(f) 5 '(f) 2 5'(f)/'(f) 

/(r)s(f) 2 + s 2 f(f)g(f) ~ 2ro(f) " 45(f) 2 + 45(f)/(f) 
5'(f) 5 '(f) 36'(f) 5 '(f) 3 ff '(f) 2 5"(f) 3r 2 i(°' 1 )(f,z)C , (°' 1 )(f,2) 3r 2 C(°' 1 )(f , z) 2 



(20) 



45(f)g(f) 46(f) 5 (f) 4 5 (f) 2 25(f) L 2 a(r )/(f )g(f) L 2 a(f)f(f)g(f) 

r 2 BW(r,z) 3g'{f)A^°\f,z) B^(r,z) a'{f)B^(f,z) 
L 2 a(f)f{f)g(f) 2glf) r + 25(f) 

f'(?)B^°\?,z) g'{?)B^\?,z) W(f)C^°\f,z) 3g'{f)C^°Hf,z) 

27(f) W 26(f) <?(f) 

3i( 1 ' 0) (f, Z )C( 1 '°)(f,z) -a^Cf.^ + B^Cf.z) = 0, (21) 
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_ 2e 2B(f,z)-2C(r,z) r 2 5e 2B(f,z) r 2 g , ^ a' {f)g' {f) b'{f)g'(r) 



f(f)g(f) 2 s 2 f(f)g(f) 2rg(r) Aa(f)g(r) 4b(r)g(r) 

f'(r)g'(f) g'(f) 2 g"(f) r 2 ^ - 1 '^, z)^ ' 1 )^, z) 3r 2 C^°^(r, z) 



4f(f)g(r) 2g(rf 2g(f) L 2 ~a{f)f{f)g{?) L 2 d(r)f(r)g(r) 

r 2 &°V{?,z) g'(f)A^{r,z) C^°>(f,z) h'{f)& lfi \f,z) b'{f)& lfi \f,z) f'(f)C^(f,z) 
L 2 ti(?)f(?)g(?) ' 2g(f) r ' ' 2a(f) ' ' 26(f) " " 2/(f) 

+ WWCWHr, z) + i(1>0) + 3(5(li0) 2 + ~ m 

2g{r) 

The index (n, m) above the functions A, B, C refer to the order of the derivative in f and z respectively while 
the primes over the functions /, a, b denotes the derivative with respect to f. The function g(f) is defined as 
g(f) = r 2 i +f 2 . 

We integrated the system of non-linear partial differential equations for some values of the cosmological 
constant with the solver Fidisol [53] based on a Newton-Raphson method, supplemented with the following 
boundary conditions 

d f A(0, z)=0, dfC(0, z) = 0, 5(0, z) - 1(0, z) = d a 

A(oa,z) = , B(oo,z) = , C*(oo, z) = 0, (23) 
d z A(f, 0) = , d z B{f, 0) = , d z C(r, 0) = 0, 
d z A{f, L) = , d z B(f, L) = , 9 z C*(f , L) = 0, 

L being the critical length given by the first order analysis. L — lix /k. These boundary conditions arc similar 
to the one used in the asymptotically locally flat non uniform black string problem [24j [25] . 

The parameter do is related to the temperature of the non uniform black hole T^ u according to 

T% u = e' do T^, (24) 

Tjj being the temperature of the background uniform solution. 

Note that here we keep the length L fixed and vary the temperature, leading to a non uniform solution 
with the critical length and different temperature than the background uniform solution. Note also that the 
entropy S NU changes along the non uniform branch and is given by 

S NU =S U f L e S ^ +d ^dz, (25) 



S u being the entropy of the background uniform phase. 

By construction, equations ([2"0"]) . ([2~1]) and ([2"2"]) admit the trivial solution A = B = C = for do = 0. This 
was checked numerically providing a crosscheck of our equations. We started with the trivial solution as an 
initial guess and gradually increased the value of do, but the solver failed to provide a convincing solution. 
Instead, starting with a combination of the trivial solution and of the first order solution and a small but non 
vanishing value of the parameter do leads to a consistent results, which are numerically robust. It has been 
checked that in the limit where do goes to zero, the trivial solution is recovered, connecting the non-uniform 
black string phase to the uniform black string. When do is increased, the solutions develop a more and more 
pronounced extremum for z = L/2, suggesting a transition to a localised black hole; this is shown in figure 
[3J where we plot an embedding of the horizon in R 3 for A = —0, 1 and different values of do. 

We compared the value of the corrections on the horizon quantities coming from the perturbativc approach 
and the same values computed from the non perturbativc approach and it appeared that our results are 
consistent, providing a crosscheck of both the perturbative and non perturbativc solutions. This is shown in 
figure [4] for three values of the cosmological constant. 

We computed the deformation parameter [11] A = Rmax/Rmin ~ 1 5 with R ma x = max{rhe ct " r '" z ^) , 
Rmin = min(rhe c ( rh ' z ' > ). It appeared that the value of A increases quicker with do for larger negative values 
of the cosmological constant; this is shown in figure [SJ 



9 



d„=0,04 




d n =0,07 






Figure 3: Embedding of the horizon surface in a three dimensional Euclidean space for ru = 1, A = — 0, 1. 
Each points on the surface represents a S3 sphere. The equation of the surface in cylindrical coordinate 
(p, 9, z) is given by p = rhe c ^ rh ' z ^ L \ 



( 5T NU /SS^ U )(SS U /ST U ) 



A=-0,1 



A =-0.2 



0.1 
a 



Figure 4: The ratio of the specific heat (5S/STh) computed with the perturbative approach and with the 
non perturbative approach. In the figure, NU refers to the non perturbative solution while U refers to the 
perturbative solution. For small values of do, the ration is close to one. 
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0.04 0.06 0.08 0.1 



Figure 5: The value of A as a function of the parameter do for = 1, L = 2tt /k, where k is the critical 
wavenumber and different values of A. The value of A increases quicker with do for larger negative values of 
the cosmological constant. 

It must be stressed that the numerical investigation was plagued by several difficulties, in particular, 
the solver was very sensitive to the mesh and failed to provide reasonable solutions once the value of the 
cosmological constant is large, this can be due to the fact that the value of the background solution becomes 
very high for large values of the cosmological constant, forcing the program to deal with extremely small 
numbers for the non uniform fields (recall that they go to zero asymptotically) and extremely large numbers 
for the background fields (which diverge like r 2 ). This is the reason why we could even with many difficulties 
investigate the solution up to only A = —0, 3. Note also that the precision of the solutions is of the order of 
5% with a 40 x 81 mesh for z x f, uniform in the z direction; it is possible to reduce this error by considering 
a denser mesh. However, it has been checked that the results presented here are compatible with results 
obtained from a denser mesh (within the numerical error) . Let us finally mention the fact that the mass and 
tension can be in principle computed since the asymptotic of the full solution should follow the Fefferman 
Graham expansion. However, it turned out that the useful subleading coefficients were impossible to extract 
in our case. However, it should be possible to integrate the Smarr relation and the first law, but we didn't 
have data enough in order to extract precise values of the mass and tension. 

8 Conclusion 

We have integrated the non uniform black strings equations in AdS up to second order in perturbation theory 
for many number of dimensions. We investigated the thermodynamical properties of these perturbative non 
uniform solutions and found evidences for the existence of a new stable phase of non uniform black strings. 
This new phase has the property that the non uniform black string length is of the order of the AdS radius. 
We propose the terminology long and short AdS black strings in contrast with the small and big AdS black 
strings terminology. The small (resp. large) black strings are characterised by a small horizon radius - AdS 
length ratio (resp. of order one and larger); the short (resp long) non uniform black strings are characterised 
by a small length in the extra direction - AdS radius ratio (resp. of order one and larger). This new phase 
is not present in the uniform case where all the thermodynamical quantities can be defined per unit length, 
so that the length plays a spectator role. For non uniform black strings, the length of the extra direction 
appears in a non trivial way in the definition of the thermodynamical quantities. 

We were not able to investigate the problem of (^-depending) critical dimension here (see [26] for asymp- 
totically locally flat spacetime). The problem comes from the extraction of numerical factors needed in 
order to construct the conserved global quantities. A possible alternative was to construct these quantities 
by integrating the first law of thermodynamics; this works qualitatively well except in the region of small 
horizon radius- AdS length ratio. Unfortunately, this is precisely the region needed. A reconsideration of the 
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numerical technique is needed in order to solve this problem. 

However, we believe these new phases of non uniform black string to be strongly related to the occurrence 
of critical dimensions in canonical ensemble. Moreover, the change of sign in the temperature correction 
should translate in a change of sign in rescaled mass correction in a rescaled mass-relative tension diagram; 
resulting if true in a connection with critical dimension in microcanonical ensemble. 

Another interesting aspect of these new phases is the connection with the boundary conformal field 
theory. Our results suggests the existence of a stable non uniform long configuration in a S 1 x S d ~ 3 spatial 
background in the CFT side. 

We also gave some preliminary results on the full nonlinear problem's solution. Despite of our efforts, 
we faced many difficulties in the investigation of the solutions, especially for larger negative values of the 
cosmological constant. However, our non pcrturbative results confirm the perturbative solutions, and indicate 
a possible transition to a localised black hole system. Unfortunately, we could not check the existence of 
the thermodynamically stable phase, predicted by the perturbative analysis. Clearly, more efforts should be 
done and a systematic investigation of the set of parameter should be investigated in the non perturbative 
regime. 

Finally, let us give some interesting perspectives to this work. One is the hunt of critical dimensions. 
An interesting question to address would be to investigate whether these sable long non uniform black 
strings could be the endpoint of the unstable AdS black string with corresponding length decay This should 
be numerically easier than problems involving topology changes and might be of peculiar relevance in the 
context of AdS /CFT correspondence. 
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